Quantum phase transitions in the Kondo-necklace model 
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Kondo-necklace model can describe the magnetic low-energy limit of strongly correlated heavy fermion ma- 
terials. There exists multiple energy scales in this model, each of which indicates a phase. Here, we study 
quantum phase transitions between these different phases, and show the effect of anisotropies in terms of quan- 
tum information properties and the vanishing of energy gap. We employ the perturbative unitary transformations 
to calculate the energy gap and spin-spin correlations for the model in spatial dimensions d — 1, 2, and 3 as well 
as for the spin ladders. In particular, we show that the method, although being perturbative, can detect the phase 
transition point by imposing the spontanous symmetry breaking, in good agreement with numerical and Green's 
function methods . We also use concurrence, a bipartite entanglement measure, to probe the criticality of the 
model. This is based on finding points or regions in which the derivative of concurrence shows non-analyticity. 
This concurrence is also able to detect critical points in indicative of that the quantum correlations in this model 
are basically bipartite. Moreover, we discuss the crossover from the 1- to 2-dimensional model via the ladder 
structures. 
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I. INTRODUCTION 

Strongly correlated systems have many competing energy 
scales driving system into a variety of phases. At zero temper- 
ature two different phases are separated by a quantum critical 
point, where the quantum fluctuations become strong leading 
to a quantum phases transition (QPT),' where a macroscopic 
change occurs in the system. The role of correlations, either 
classical or quantum, between underlying degrees of freedom 
is the key to understanding QPT. The latter is called entan- 
glement, which is a measure of nonlocal correlations of quan- 
tum manybody states. Entanglement is a versatile quantity in- 
tensively studied in quantum information theory and strongly 
correlated systems providing an exchange of ideas between 
two apparently distinct fields. Indeed, the entanglement has 
been recognized as an essential resource in quantum compu- 
tation and communication 2 . Entanglement based ideas lead 
to distinguish criticality from off-criticality in ld», area law 
scaling of entanglement entropy 4 and topological sub-leading 
term 5 ^. 

Kondo lattice modeli is one of the promising models in 
the study of strongly correlated heavy fermion compounds 8 , 
where the partially filled shells (f) are screened by outer shells 
(s or p). Inner shells are localized and can be effectively char- 
acterized by magnetic impurities located in the electron gas of 
outer shells. In the presence of strong interaction between the 
magnetic impurities and electron gas, a more simplified ver- 
sion of Kondo lattice, namely Kondo-necklace (KN) model 9 
arises (see SecHH]). The phases of the Kondo-necklace model 
in spatial dimensions d=l,2,3 have been studied by different 
methods. While numerical Monte Carlo^ and Density Ma- 
trix Renormalization Group' 1 (DMRG) results confirm the ab- 
sence of QPT in Id, a phase transition between antiferromag- 
netic (AF) and Kondo singlet phases occurs in higher dimen- 
sions. 

In this work we study QPT of the anisotropic Kondo- 
necklace model in spatial dimensions d = 1, 2, and 3 by 
employing perturbative continuous unitary transformations 



(PCUT)! 2 .. As the extreme limit of the spin liquid phase of the 
model is characterized as isolated dimers, the model is well 
suited through the above method, then the inter-site interac- 
tion is treated perturbatively. Moreover, the system is studied 
at the thermodynamic limit. Other merits of the method are as 
follows. At the first glance, with implementing PCUT on the 
contrary of some numerical methods, we can obtain energy of 
states analytically as a function of Hamiltonian parameters. 
It helps us to find some other underlying phase transition wit- 
nesses such as concurrence (from correlation functions that 
they themselves can be calculated from derivatives of ground 
state energy correspond to Hamiltonian parameters at zero 
temperature). In addition, in this method we don't need know- 
ing eigenstates of the system for finding ground and excited 
states energies and we can calculate energies independent of 
finding systems eigenstates (In some numerical methods such 
as quantum Monte Carlo we should guess a prime trial state). 
Last, in this method, we don't have any size and dimension 
restrictions that put lots of obstacles in front of some other 
numerical methods. We can use this method for arbitrary di- 
mensions and system sizes. 

We further characterize the QPT by using "concurrence"— 
as a measure of entanglement between a pair of spin- 1/2 par- 
ticles. We will use entanglement of the ground state in order 
to characterize the phase transition. The entanglement prop- 
erties of the one dimensional Kondo-necklace model with a 
few sites in order to have analytical results has already been 
addressed-^. In the current work we will take the thermo- 
dynamic limit via PCUT and probe the phase transitions in 
spatial dimensions d = 1, 2, and 3. Despite its local nature, 
nonanalytic behavior of derivatives of concurrence signals a 
phase transition in the system, and its scaling in the vicinity 
of critical points is connected to the universality class of the 
mode l 1516 . We will show that the critical points detected by 
such nonanalyticities are not only in well agreement with the 
gap closing points but with mean-field theory^— and Green 
function 2 ^ results. Moreover, we discuss the crossover from 
the non-critical one-dimensional model to the critical two di- 
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mensional one investigating the properties of the model on the 
ladder geometry. 

The paper is organized as follows. In Sec. [TT] the PCUT 
method is briefly reviewed. Next, in Sec. [Ill] the explicit ex- 
pression for ground state energy and excitation spectrum of 
the Kondo-necklace model are derived by PCUT method, and 
the QPT of the model is discussed in Sec. [TV] Finally, we 
present the summary and discussion of our results. 



for Hq. Then, the generator of transformation is obtained or- 
der by order in terms of perturbation parameter A. Moreover, 
the off-diagonal part evolves to reach the I = oo effective 
Hamiltonian. However, PCUT has two prerequisites: 

1. The unperturbed Hamiltonian (Ho) has to be composed 
of equidistant energy levels which are bounded from be- 
low. 



II. PERTURBATIVE CONTINUOUS UNITARY 
TRANSFORMATION 

There are few examples in strongly-correlated systems 
where one can find exact solutions, i.e., exact ground state 
and excitation spectrum in addition to some correlation func- 
tions. However, numerical methods give an approximation for 
the physical properties of these systems. One of these meth- 
ods, which can diagonalize a Hamiltonian properly, regardless 
of system size, is the continuous unitary transformation. This 
method had been introduced simultaneously in Refs. I2H - I23I 
for condensed matter systems and quantum chromodynamic 
problems. In this approach the Hamiltonian is considered as 
a function of a flow parameter (0- The Hamiltonian is trans- 
formed to a simpler form (diagonal or band-block diagonal) 
under a flow equation which is based on applying infinite 
numbers of infinitesimal unitary transformation. 

Let us define the Hamiltonian, H{1), as a function of a con- 
tinuous parameter I where H(0) is the bare Hamiltonian and 
i? (00) is the final diagonalized effective one. The evolution 
of the Hamiltonian in terms of the continuous unitary trans- 
formation is given by the following sets of equations 



2. The perturbing Hamiltonian (H s ) should be written as 



H(l) = U{1)H(0)U{1)\ 
dH{l) =[,«), H(l)], 



dl 



(1) 



transformation and the anti- 



— W is the generator of this 



where U(l) is the unitary 
hermitian operator r)(l) = 
transformation. There are several choices for rj(l) which de- 
pends on the Hamiltonian and the possibility to find a closed 
form for the transformation. For instance, Wegner considered 
H(l) = Hd{l) + H n d{l) where Hd is the diagonal part of 
Hamiltonian in a specified base ket and H n d is the correspond- 
ing off-diagonal part. Accordingly, the generator is defined by 



ri(!) = [H d {0),H(!)] 



(2) 



However, it does not lead to an analytic solution for I —> 00, as 
a general case. Moreover, the Wegner choice mixes the blocks 
of Hamiltonian during transformation if the original Hamilto- 
nian is block diagonal. To overcome this problem Mielke 24 
proposed a generator to preserve the band diagonal form of 
the initial Hamiltonian and accordingly Knetter and Uhrig in- 
troduced the perturbative continuous unitary transformation 
(PCUT)ii. In this approach the original Hamiltonian is de- 
composed to two parts H (0) = Hq + XH S , where the unper- 
turbed Hamiltonian is Hq = Hd(0) and H s represents the off- 
diagonal part at I = which is considered as a perturbation 
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T n for finite value of N, where T n in- 



creases or decreases the energy quanta of the specified 
levels defined by Hq, i.e., [Hq, T n ] = nT n . 

We follow the formalism defined by Uhrig and Knetter, in 
which the generator of transformation is defined by 



V (l) = J2* k E ^n{M{m))F{km)T{m), 

k—1 \m\— k 
k 

M(m) = J2 m i' (3) 

i=l 

where k is the order of perturbation, F(l;m) are real val- 
ued functions which are obtained by recursive differential 
equations^ and m denotes a sequence of quanta labels 

m = (mi,m 2 ,m 3 ,...,m k ), 
m, G {0,±1,±2,..-±JV}, 



T ( tti ) — T m 1 T m 2 T m g . . . T mfc . 



(4) 



Moreover, |m| = k means the summation is done on all pos- 
sible configurations of a sequence with k members. In the 
above equation, M(m) is the total number of energy quanta 
which is created or annihilated by T(m). In this respect, the 
Hamiltonian is given in the following form 



Hi=H + J2^ k Yl F ^™) T (™)- 



(5) 



ml— k 



In order to get a final diagonalized Hamiltonian, the num- 
ber of energy quanta which is created or annihilated at I = 00 
must be zero. Therefore, we can write the effective Hamilto- 
nian as 



C(m)T(m), (6) 



H eff =Ho + J2* k E 

fc=l \m\=k,M(m)=0 



where C(m) = F(l = 00; m). The solution of the flow equa- 
tion (Q]) for the generator given by Eq. (0 and the Hamiltonian 
defined in Eq. (0 leads to the coefficients C(m). Hence, the 
effective Hamiltonian is obtained perturbatively as a sum of 
different orders of perturbation. 
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Moreover, the spectrum of a local dimer is composed of a 
singlet and triplet state representing two quanta, and 1, re- 
spectively. Thus the effect of H s on two neighboring singlets 
changes the quanta by n with n = 0, 1, 2 which provides the 
same coefficients (C(m)) derived in Ref. [l2| for our calcula- 
tions. 



T. 



'i+1 



FIG. 1. One dimensional Kondo-necklace model. Solid and grey 
circles stand for spins of itinerant and localized electrons. 



III. PCUT FOR KONDO-NECKLACE MODEL 

The anisotropic Kondo-necklace model on a hypercubic lat- 
tice is defined by the following Hamiltonian, 



y y 



At 



(id) 



(7) 



where Tj is the spin of itinerant electrons on the d-dimensional 
hypercubic lattice with L d sites, Si is the spin of localized 
electrons connected to each lattice point and represents 
the nearest neighbor sites on the lattice, and Jo and A are, 
respectively, the local exchange coupling and anisotropy pa- 
rameter. Here J and -q^ are the exchange and anisotropy cou- 
plings on the lattice, respectively. The one-dimensional lattice 
is shown in Fig.Q]and a generalization to two and three dimen- 
sional system can be deduced from this figure. We have con- 
sidered the fully anisotropic case, which enables us to study 
the anisotropy effects as well to obtain the correlation func- 
tions via introducing a generating function. For simplicity, we 
define J/ Jo = A and set Jo = 1. To apply PCUT we consider 
the following decomposition of Hamiltonian 



H — Hq + \H S . 
H = E s » ■ T i} 



(8) 



<id> 



where a = (A — 1)/A. The anisotropic part of local inter- 
actions are considered in the perturbing part (H s ) in order to 
have an equidistant spectrum for Ho which satisfies the first 
criterion for PCUT approach. The ground state of H is a 
Kondo-singlet state for Jo ^ J which justifies A = J/ Jo 
to be small. In other words, we trace the quantum phase tran- 
sition from the Kondo-singlet phase where A is a small pa- 
rameter. The Kondo-singlet state is the ground state of the 
unperturbed Hamiltonian (Ho) which is a product of singlets 
on local dimers (r^-Si pair). 



A. T„ operators 

To find Heft we need to calculate the T(m) operators which 
essentially define the effect of H s on a base ket. As far as 
the Hamiltonian is composed of two body interactions the ef- 
fect of H s can be represented by its effect on two neighboring 
dimers, i.e., 



He 



E*i 

(id) 



(9) 



A dimer is a pair of Tj spin and its corresponding neighbor- 
ing impurity spin Sj which can be in a singlet or triplet state 
configurations. Accordingly, the effect of if J J on two neigh- 
boring dimers may change the quanta of the state by 
n = 0, ±1,±2 which is represented by q n . Therefore, the 
pair Hamiltonian of H s is written in the following form: 



= q- 2 + q-i +qo + qi + 



?2- 



(10) 



Two neighboring dimers are represented by 16 configurations 
depending on the configuration of each dimer to be in a sin- 
glet or triplet state. We have summarized the effect of all q n 
operators on all possible 16 states of two neighboring dimers 
in Table |A]of AppendixlAl. The translational symmetry of the 
Hamiltonian allows us to use TablelAlfor any pair of neighbor- 
ing dimers. According to Eq. (O the perturbing Hamiltonian, 
H s is written in terms of a sum on q n operators. For instance, 
in the one-dimensional model which has L dimers with peri- 
odic boundary condition we arrive at 



i— 1 i—l i—1 i=X i—1 



Ti+T 2 . 



; = 1 
(ID 



In the above equation, T n is the sum of q n operators on all 
bonds which connects two neighboring dimers on the lattice. 



B. Ground-state energy 

According to the PCUT approach the ground-state en- 
ergy (GSE) is given by (0|Jf e ff|0), where H e s is given in 
Eq. ©, and |0) is the direct product of singlet states over 
all dimers. We have calculated the ground-state energy per 
spin (e = Eq/ (2L d ) for the fully anisotropic model and on 
e?-dimensional hypercubic lattice up to the 4th order of per- 
turbation. Note that there are two spins (r, s) corresponding 



4 



TABLE I. The ground-state energy per site for the one-dimensional 
KN model with r] x = rj y = l,rj z =0 and A = 1. A comparison 
between the 4th order PCUT and Lanczos results is given by the 
relative error |eLanczos — epcurl/eL 



A 


^Lanczos 


epcuT 


relative error 





-0.37500000 


-0.37500000 





0.1 


-0.37531250 


-0.37531250 





0.2 


-0.37624999 


-0.37625000 


0.00000003 


0.3 


-0.37777517 


-0.37781250 


0.0001 


0.4 


-0.37999946 


-0.38000000 


0.000001 


0.5 


-0.38281007 


-0.38281250 


0.000006 


0.6 


-0.38631563 


-0.38625000 


0.0002 


0.7 


-0.39040284 


-0.39031250 


0.0002 


0.8 


-0.39496177 


-0.39500000 


0.0001 


0.9 


-0.40020219 


-0.40031250 


0.0003 


1 


-0.40611293 


-0.40625000 


0.0003 



with the Lanczos ones where the maximum error is less than 
0.1% for A = 1. It confirms that the PCUT results are reliable 
as far as A < 1. We have also made a similar comparison 
for other values of anisotropies which lead to the same con- 
clusion. We have also calculated the ground state energy of 
Kondo-necklace model on ladder geometry which is a quasi 
one dimensional model. A three legged ladder is shown in 
Fig. |2] This is important to study the crossover between one 
and two dimensional systems. The n-leg KN ladder is a set 
of n KN chains which interact via the r-spins as given by the 



1=1 



- T S s s s 



1=3 







to each lattice point, r on the hypercubic lattice and s a lo- 
cal spin moment. Moreover, the following result is valid for a 
lattice in the thermodynamic limit (L oo): 



FIG. 2. Three legged Kondo-necklace ladder. Solid and grey circles 
stand for spins of itinerant and localized electrons. 



(2ad)- 



dX 



2L d 

+ ^\ 3 (-3r lx r ]y T lz + 2(2ad){r)l + rfi)) 

+±t((4d - iM + r,l + a 2 - iv 2 M + v 2 y ) 

-4r) 2 x r,*-4(4d-3)(ri$ + r)* + 4) 

-I6(2ad) 2 {?i + rfi) + 32(2ad)r lx r ]yVz ) , (12) 



where A = A/4. The first term in the above equation is the 
singlet ground state energy, i.e. the zeroth order value. We 
have compared the PCUT results with the Lanczos diagonal- 
ization ones for the one-dimensional isotropic case in Table Q] 
The Lanczos computations have been done on several lengths 
up to 24 spins (L = 6, 8, 10, 12) where we found a tiny finite 
size dependence which concludes that the 24 spins model rep- 
resents the thermodynamic value. The last column in Table 
U shows the percentage error of the PCUT results compared 



following Hamiltonian 

L n 

i=l 3=1 
L+l n 



(13) 



i=l j = l 
L n-1 

J' ( / ;( r lx T i3' r f,.i+l + T ly T ij T i,j+l + Vz T ij T tj+l))i 



i=l 3 = 1 

where the local term is isotropic, the anisotropy parameters 
along the chain direction is given by r\ a and the anisotropy 
along the rung direction is represented by rf a . The exchange 
interaction along the rungs is J' and its ratio to the leg inter- 
actions is given by x = J' / J. The periodic boundary condi- 
tion is assumed along the leg direction while it is open along 
rungs. The ground state energy per spin for n-leg KN ladder 
to the 4th order of perturbation expansion in PCUT approach 
is obtained to be 



_ladder_ ^0 



ladder 



+ 



A 4 



_ X_ 

8 ~ T 
(n-1) 



1G 

vi - % - Vz 



I 2 i 2 i 2 \ i 2 ( n 1) / '2 
Wx + Vy + Vz) + X 

n 

L; (vl + v 2 v + v 2 z)(v' x 2 + v y 2 + Vz 2 ) + ^ [ - 



y 2 + €) 
- i 



3A_ 



VxVyVz 



■VxVyV, 



(fix + V + Vz) + HVxVy) + KVxVz) + HVyVzY 



4 ~ 4 - 4 i ^(32 X 2 +^)(^ + n ^ + ^)-^ i -^^^^^ 



(14) 
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K. 

FIG. 3. Contour plot of gap function G2(k x , k y ) of 2d KN model. 
The gap gets minimum at (tt, 7t). The corresponding parameters are 
A = 0.5, rj x = rj y — 1, r\ z = 0, and a = 0. 



The above equation reduces to the energy of d = 1 dimen- 
sional KN model for n = 1 and also for m = 0. It should be 
noted that Eq.dl~4b has been calculated for the isotropic local 
interaction a = 0. 



In contrast to the ground state energy the spatial dimension 
does not enter directly to the dispersion of magnons, because 
the calculation depends on the position of the triplet excited 
dimer. A similar calculation leads to the magnon dispersion 
for two and three dimensional hypercubic lattice which is pre- 
sented in AppendixlBl 



The minimum of the magnon energy defines the energy gap 
which appears at the antiferromagnetic wave vector k x = it 
for the one dimensional model and gives the following depen- 



C. Excitation spectrum 

The one magnon dispersion can be obtained using PCUT 
approach. A magnon is created by exciting a singlet dimer to 
a triplet in the Kondo-singlet ground state. We show the one 
magnon state by \j) =| s, s, • • • ,tj,- ■ ■ ,s) where the singlet 
dimer at position j has been replaced by a triplet state. Among 
the three states in the triplet set (t ± , t°), the energy of is 
lower than t°. Thus we choose either t + or t~ for constructing 
the one magnon state. The Hamiltonian conserves the number 
of triplets, so all states with a single triplet have the same en- 
ergy. Therefore, the linear combination of all one magnon 
states at different positions is defined as the elementary ex- 
cited state which is represented by |k) for a d-dimensional 
hypercubic periodic lattice, 

|k) = -^=^e lk - r | SS 8- -88). (15) 

The dispersion of the magnon spectrum is calculated by the 
following equation 

w(k) = (k|ifeff|k) - E Q . (16) 

Due to the existence of a triplet dimer in the excited state, 
the amount of calculation of Eq. dT6l increases dramatically 
for the 4th order of perturbation. Therefore, we had to re- 
strict our calculation to the 3rd order of perturbation for the 

one magnon spectrum. The magnon dispersion of the one di- 
mensional model is expressed by the following relation with 

T}x = Vy = 1. 



(17) 

I 

dence of the gap on the model parameters, 

Gi = \ + j + A[-2] + A 2 [| - A A + 4 Vz + ^] 

+ A 3 [4+^ + ^ 2 ]. (18) 

The gap, G\, is always nonzero for the interested range of 
parameters representing no magnetic order for the one di- 
mensional model while G*2 and G3 behave differently. The 
magnon dispersion of the two dimensional lattice is plotted in 
Fig. ([3]) for A — 0.5, r\ x = r) y = 1, f] z = and a = 0. The 
excitation energy gets its minimum at k = (ir,ir) which ver- 
ifies the antiferromagnetic ordering. For the two dimensional 



wi(k) = 1 + A[4a + 2cos(fc T )] - A [r\ 7 2 - 477 z cos(fca;) + cos(2fc x )] 



T 3 r 17 
+ A [--a 



13 
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+ 4acos(2fc. E ) — 5cos(k x ) + cos(3fc K ) + — -r\ z + 2r] z cos(2fca;) — —r\ z cos^)]. 



lattice, the energy gap is given by 

G, = A-4A-A 2 [8?7 2 + 2^ + 10A-4] + A 3 [21?72+40r/ z -l]. 

(19) 

A similar calculations show that for three dimensional lat- 
tice the excitation energy is minimized at the antiferromag- 
netic vector k = (it, n, tt) where the energy gap is 

G 3 = -\ + - 6A + A 2 (-| - 12„. + 3rf z + |a) 

+ A 3 (ip-69 % + H\ 2 ). (20) 

The analysis of G2 and G3 are given in the next section where 
quantum critical point (QCP) is recognized at the position of 
zero gap. 




k 

FIG. 4. Energy gap (G2) vs. A for different values of A for rj x = 
rj y = 1 and r\ z = on the two dimensional lattice. 



IV. QUANTUM PHASE TRANSITION 



The quantum phase transition in KN model is a competi- 
tion between the Kondo-singlet phase and the antiferromag- 
netic long range order. The U(\) symmetric one dimensional 
model is always in the Kondo-singlet phase for any value of 
A 17 ' 25 which agrees with the nonzero gap (Gi) obtained in 
Eq. ( fT~8b - However, the one dimensional model with Z2 sym- 
metry shows a quantum phase transition to antiferromagnetic 
order at A c = 2.22 for the Ising interaction between r-spins, 
i.e., -q x = rjy = and r\ z = 1—. Unlike the one dimen- 
sional model, both two and three dimensional KN Hamilto- 
nian with U(l) symmetry show a quantum phase transition 
from the disordered Kondo-singlet to the antiferromagnetic 
ordered phas o 17 ' 20,25 . The quantum critical point depends on 
the anisotropy parameters; however, the critical value of the 
control parameter is less than one (A c < 1) for both 2- and 
3-dimensional models. It supports our approach to study the 
quantum phase transition of KN model by PCUT 

We study the QPT of KN model by implementing two dif- 
ferent measures. (A) The energy gap which vanishes at the 
QCP. (B) The derivative of bipartite entanglement that shows 
a singular behavior at the QCP. It should be emphasized that 
PCUT gives a results for the thermodynamic limit (L — > 00) 
although its accuracy is given by the order of perturbation. 



A. Energy gap 

We obtain the quantum critical point (QCP) of the two and 
three dimensional KN model for different anisotropy parame- 
ters. We have plotted the energy gap of the two dimensional 
model, Eq. ( fT9b , in Fig. [4] and the three dimensional gap, Eq. 
( |20t , in Fig. [5] Different plots belong to different anisotropy 
values (A) which clearly manifest the strong dependence of 
QCP on A. For the two dimensional model the results from 
other methods such as mean-filedi 7 - and Green's function 2 ^ ap- 
proaches are available, and we have summarized the variation 
of QCP on the anisotropy parameter in Table HI1 The results of 
all methods are in agreement with each other (except at small 




FIG. 5. Energy gap (G3) versus A for different values of A and r\ x 
r] v = 1, rj z = in the three dimensional lattice 



values A). It is clearly seen that as the anisotropy parameter 
A on the bonds increases the phase transition at critical point 
A c occurs at higher values. This is sensible because the gap 
between singlet and triplet states increases with A. Thus a 
higher value of A c is needed to close the gap and drive the 
Kondo singlet phase to the antiferromagnetic phase. 



B. Concurrence 

In this subsection the entanglement between spin degrees of 
freedom is used to detect the phase transition point. The role 
of entanglement in characterizing QPT has been studied in the 
past few years. Reference j_15f showed that the concurrence as 
a measure of entanglement, though being a local quantity, can 



TABLE II. The quantum critical point for different anisotropy (A) 
in two dimensional model. PCUT results are compared with the 
corresponding from mean-field (MF) and Green's function (GF) ap- 
proaches. 



A 


0.2 


0.4 


0.6 


0.8 


1.0 


A c (PCUT) 


0.2 


0.40 


0.56 


0.68 


0.77 


A c (MF) 


0.42 


0.50 


0.56 


0.63 


0.70 


A c (GF) 


0.39 


0.46 


0.52 


0.58 


0.65 



be considered to trace the QPT. They found that the deriva- 
tive of concurrence with respect to the control parameter di- 
verges at the QCP for the Ising model in a transverse mag- 
netic field. Finite size scaling and critical exponents coincide 
with the universality class of the phase transition. The non- 
analytic behavior of quantum correlations which are encoded 
in entanglement has been discussed in several works for spin 
models2£ii. as well as itinerant system s 3 ^ - —. In this respect 
we study the derivative of concurrence for the KN model to 
find the QCP and its dependence on the anisotropy parame- 
ters. In particular, we use concurrence as a profound measure 
of entanglement between two spin 1/2. 13 

The concurrence of a pair of spin- 1/2 is obtained by the 
following expression, 



(21) 



C = max{0, v Ai 



A 4 }, 



where Ai < A2 < A3 < A4 are the eigenvalues of the matrix 
PabPab, in which pab is the reduced density matrix for the 
pair labeled by A and B, 



PAB 



y * v 
' °bPab°a 



(22) 



The reduced density matrix can be written in terms of two 
point correlation functions, 



PAB 



°A< 



B • 



P, 



'Bh 



(23) 



where i, j runs over all Pauli matrices and unit operator in the 
corresponding Hilbert space. The \<Ja{b) can be either s or 
t spins defined in the KN Hamiltonian. We define the follow- 
ing generating functions to calculate the correlation functions 
from the ground state energy obtained in previous sections, 



T; S , 



1 



d 



d\L d drip 

_L d 

~ d\L d d~c 



H. 



-H, 



(24) 



and taking into account the normalization of the ground state 
we will arrive at 



KTiT- 



1 



d 



d\L d drip 

_L d_ 

~ d\L d da 



En 



En 



(25) 




FIG. 6. Concurrence of two neighboring spins of conduction elec- 
trons (n, Ti+i) for the 1, 2, 3d isotropic KN model and n-leg ladder 
versus A for r) m = 1, r] z = 1, rj y = and A = 1. The ground state 
preserves Z2 symmetry which is equivalent to no symmetry break- 
ing. 



The explicit expressions of the correlation functions in terms 
of model parameters for two neighboring spins are presented 
in AppendixICl 

We consider two cases to study the concurrence and its 
derivative for KN model. In the first case, the ground state is 
supposed to have Z2 symmetry, i.e., no spontaneous symme- 
try breaking (SSB). We get zero magnetization for this ground 
state, m a = (r a ) = 0. Therefore, the concurrence between 
two neighboring spins of itinerant electors on the hypercubic 
lattice is given by the following equation: 



G-- 



i T i+ll 



E 

P=x,y,z 

= -^(Vx+Vy + Vz)X + ^(VxVy + VxVz + Vh)~>? 
J 3 

+_ [(44 _ i)( % + Vy + +1 f y + r,l) 

+2(VxVl + VyVl + VxV 2 z + VzVl + VyV 2 z + Vzvl) 
+4(4d-3)(4+r 1 3 y +,f z )] -i. (26) 



We have plotted the concurrence of two neighboring r-spins 
given by Eq. (l26j i in Fig. |6]for one, two and three dimensional 
KN lattice (in addition to the n-leg ladder structures which 
will be discussed latter). None of plots show a non-analytic 
behavior to signal out a quantum phase transition. We have 
also examined the derivative of plots presented in Fig. [6] and 
found no divergent behavior. Referencel35l pointed out that the 
symmetry breaking of the ground state is necessary to be taken 
into account to find the quantum phase transition in some spin 
models. Accordingly, we will calculate the concurrence in the 
presence of SSB. 

The PCUT approach is a systematic perturbative scheme to 
consider the effect of interactions on a noninteracting Hamil- 
tonian (Hq). In our model, the ground state of Hq is a Kondo- 



singlet state which preserves the Zi symmetry of the Hamilto- 
nian. And, all order perturbative results respects the sym- 
metry which gives both magnetization and staggered magneti- 
zation to be zero. The SSB can be added to our PCUT results 
through a mean-field approximation. In this case, the expec- 
tation value of a single spin is approximated from the pair 
correlation function, i.e., 

(T?T? +n )^{Tn(T? +n ) = ( m a ) 2 . (27) 

The off-diagonal correlation functions are also approximated 
as the product of the corresponding spin expectation values 
in the presence of SSB. The concurrence, Eq. (I2TI 1, is calcu- 
lated again in the presence of SSB and the corresponding plots 
for one, two and three dimensional models are presented in 
Fig. [7] The rapid change of concurrence close to QCP leads 
to a seemingly divergent behavior for the derivative of con- 
currence which is shown in Fig. [8] It justifies that the diver- 
gent behavior for derivative of concurrence which comes out 
of a ground state with SSB represents the quantum critical 
point. We have summarized the QCP obtained by this method 
in Table [TIT] along with the results obtained by the vanishing 
of energy gap and also the results of mean-field and Green's 
function approaches. The 2c? result which has been obtained 
by the derivative of concurrence is close to the the GF QCP 
while there is a difference for the result which is given by the 
zero point of energy gap. The reason is mainly due to the 
third order calculations for energy gap compared with the 4th 
order calculations for the ground state energy from which the 
concurrence has been obtained. A controversial issue rise up 
here, the U(l) symmetric one dimensional KN model shows 
a quantum phase transition as witnessed by the divergent be- 
havior of derivative of concurrence. Although a qualitative 
statement and the mean-field approac h 1719 rule out the antifer- 
romagnetic order for the one dimensional model, the present 
PCUT results predicts a QCP in agreement with the Green's 
function approach^ 6 . We will discuss more on this issue in the 
next section. 

TABLE III. The QCP of the isotropic d=2, 3 dimensional KN model 
obtained by mean-field (MF), Green's function (GF), derivative of 
concurrence in the presence of SSB by PCUT (C SSB ), zero gap by 
PCUT (G 2 = or G- 3 = 0). 





MF 


GF 


£,SSB 


Zero gap 


d = 2 


0.70 


0.65 


0.68 


0.77 


d = 3 


0.40 




0.63 


0.59 



V. CONCLUSION AND DISCUSSIONS 

We have applied the perturbative continuous unitary 
transformation to study the quantum phase transition and 
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FIG. 7. Concurrence of two neighboring spins of conduction elec- 
trons (n, Tj+i) in the presence of explicit breaking of Z2 symmetry 
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FIG. 8. Derivative of concurrence (in the presence of SSB, dC/dX) 
for the plots presented in Fig. [7] The singularities signal the phase 
transition. 

anisotropy effects on the Kondo-necklace model. The advan- 
tage of PCUT is that it gives the results for the thermody- 
namic limit (L 00) which means no finite size effect ex- 
ists in these results. However, it is a perturbative approach 
and the accuracy of results depends on the order of pertur- 
bation which has been implemented to calculate the effective 
Hamiltonian. We have obtained the effective Hamiltonian and 
consequently the ground state energy to the 4th order of per- 
turbative scheme. A systematic calculations to higher orders 
of perturbation [0(A 10 )] is possible for future investigations. 

We have calculated the ground state energy to order A 4 for 
the hypercubic KN model in d = 1,2,3 dimension and also 
for the n-legged ladder. To justify the accuracy of the PCUT 
results we have compared the ground state energy per site for 
the one dimensional model with the Lanczos results which is 
presented in TableU The Lanczos computation has been done 
on a chain up to 24 spins. The finite size scaling analysis 
on the Lanczos results shows a tiny finite size effect which 
concludes that the results for 2L = 24 represents the values 



of thermodynamic limit. The relative error presented in Table. 
Ujustifies the high accuracy of the PCUT results. 

The excitation energies of the KN model have been ob- 
tained within PCUT implementation. The dispersion relation 
for d-dimensional KN model have been obtained as an ana- 
lytic expression, Eqs. (fTTJ, dBU , and dB21 i. The minimum 
of excitations defines the energy gap which vanishes at the 
quantum critical point. The quantum critical points of the two 
and three dimensional KN model and its dependence on the 
anisotropy parameter A have been presented in Table [II] and 
Figs. [4] and [5] Both figures show that the QCP depends on 
A strongly. However, the divergent behavior of the derivative 
of concurrence seems to give more accurate results than the 
zero gap position since the concurrence contains 0(A 4 ) terms 
while the gap has 0(A 3 ) ones. 

More information on the nature of quantum phase transition 
in the KN model can be obtained by computing entanglement 
in the system. We have obtained the concurrence as a measure 
of bipartite entanglement within the PCUT approach. We have 
defined a generating function which gives the two-point corre- 
lation functions and consequently the bipartite entanglement 
of the KN model. According to Ref . [Tsl the derivative of con- 
currence is singular at the QCP. We have calculated concur- 
rence both without and with SSB. In the absence of SSB, we 
obtain zero for magnetization and also the off-diagonal corre- 
lation functions. We have observed no singular behavior in the 
derivative of concurrence versus A without SSB, which is in 
agreement with the DMRG computations reported by Ref. H3- 
This may lead us to conclude that the quantum correlations are 
beyond the bipartite entanglement. However, the derivative of 
concurrence which has been obtained in the presence of SSB 
gives a singular behavior at a specific value of A, the quantum 
critical point, see Fig. [8] It justifies the idea of Ref. [35| that 
the presence of SSB is necessary to trace the quantum critical 
points using the generalized entanglement entropy. Moreover, 
it manifests that the quantum correlations in the KN model are 
of bipartite nature. 

The calculated values of QCP for two and three dimen- 
sional models are more reliable than the one dimensional case. 
Because, A c <C 1 (A c < 1), which justifies the accuracy of 
PCUT for d — 2, 3 while it is close to one (A c < 1) for d = 1. 
Moreover, the mean field approximation fails for d = l 19 , 
which sheds doubt on using the mean field values to calcu- 
late concurrence with SSB while mean-field results have been 
reported to be reliable for d = 2, 3 cases 1 -. The mean field 
approach 17 yields a Kondo singlet phase for the whole range 
of A, i.e., no QPT while Green's function approach predict 
a QPT from Kondo singlet to antiferromagnetic order—. A 
qualitative description predicts no QPT for the U(l) symmet- 
ric d = 1 KN model because we get no long range order for 
both extreme limit of A. The A = gives the Kondo singlet 
state and A — > oo leads to the one dimensional spin- 1/2 XXZ 
model which has a spin-fluid ground state with no long range 
order. In this respect the quantum phase transition in d = 1 
KN model with U(l) symmetry is controversial. However, 
the Zi symmetric KN model shows a QPT from the Kondo 
singlet phase to antiferromagnetic ordered phased. It opens 
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FIG. 9. Ground-state energy per site for one, two, three dimensional 
hypercubic Kondo-necklace model and also the n-legged ladder. 



the question for more investigations on the one dimensional 
U(l) symmetric KN model. An interesting result which can 
be drawn from our PCUT approach is the properties of KN 
model on the n-legged ladder geometry. We have plotted the 
ground-state energy (GSE) per spin (e) for n-legged ladder 
in addition to the 1-, 2-, 3- dimensional model in Fig. (O. 
The plots of n-legged ladder are close to the value of e for 
the d = 2 dimensional KN model which hints the behavior of 
ladders is more similar to the d = 2 model than the d = 1 
case. Clearly, GSE per spin reaches the d = 2 values in the 
n — > oo limit. We observe a quantum phase transition for 
the n-legged ladder by the abrupt change of concurrence with 
SSB and accordingly the divergent behavior of its derivative 
with respect to A presented in Figs. |7]and[8] The value of A c 
for the n-legged ladder is between the corresponding values 
for d=l and d=2 cases. However, the QCP for n = 2 legged 
ladder experiences a jump from the d = 1 counterpart. It can 
be understood simply from the expressions given in Eqs. (IT~2b 
and ( TT4T > where some extra terms of the ladder GSE will be 
omitted for n = 1. Moreover, the concurrence without SSB 
shows clearly that the ?i-legged results are similar to the d = 2 
KN model which is witnessed in Fig. [6] Although the results 
of concurrence without SSB are not appropriate to extract the 
quantum critical point they are more reliable than the results 
with SSB to extract the qualitative behavior because there is 
no mean field approximation in Fig. [6] We conclude that 
the n-legged KN model has a quantum phase transition from 
the Kondo-singlet phase to the antiferromagnetic ordered the 
same as the d = 2 dimensional KN model. 



ACKNOWLEDGMENTS 

The authors would like to thank N. Ghasemi for fruitful 
discussions. A. L. would also like to thank S. Mahmoudian 
for several discussions specially on the ladder structure. This 
work was supported in part by the Center of Excellence in 



10 



Complex Systems and Condensed Matter (www.cscm.ir). 



* langari@sharif.edu, http://spin.cscm.ir 

1 S. Sachdev, quantum phase transition (Cambridge University 
Press, Cambridge, 2000). 

2 M. A. Nielsen and I. L. Chuang, Quantum Computation and 
Quantum Communication (Cambridge University Press, Cam- 
bridge, 2000). 

3 P. Calabrese and J. Cardy, J. Stat. Mech. P06002 (2004). ibid, J. 
Phys. A: Math. Theor. 42, 504005 (2009). 

4 J. Eisert, M. Cramer, and M. B. Plenio, Rev. Mod. Phys. 82, 277 
(2010). 

5 M. Levin and X.-G. Wen, Phys. Rev. Lett. 96, 110405 (2006). 

6 A. Kitaev and J. Preskill, Phys. Rev. Lett. 96, 1 10404 (2006). 

7 H. Tsunetsugu, M. Sigrist and K. Ueda, Rev. Mod. Phys. 69, 809 
(1997). 

8 A.C.Hewson, The Kondo problem to Heavy Fermions (New York, 
Cambridge University Press, 1993). 

9 S. Doniach, Physica91B, 231 (1977). 

10 R. T. Scalettar, D. J. Scalapino, and R. L. Sugar, Phys. Rev. B 31, 
7316 (1985) and references therein. 

1 S. Moukouri, L. G. Caron, C. Bourbonnais, and L. Hubert, Phys. 
Rev. B 51, 15920(1995). 

12 C. Knetter and G.S. Uhrig, Eur. Phys. J. B 13, 209 (2000). 

13 W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). C. H. Bennett, 
D. P. DiVincenzo, J. A. Smolin, and W. K. Wooters, Phys. Rev. A, 
54, 3824 (1996). 

14 A. Saguia, M.S. Sarandy, Phys. Rev. A 67, 012315 (2003). 

15 A. Osterloh, Luigi Amico, G. Falci and Rosario Fazio, Nature 
416, 608 (2002). 

16 M. Kargarian, R. Jafari and A. Langari, Phys. Rev. A 76, 60304 
(R) (2007). ibid, Phys. Rev. A 77, 032346 (2008) 

17 A. Langari and P. Thalmeier, Phys. Rev. B 74, 024431 (2006). 

18 P. Thalmeier and A. Langari, Phys. Rev. B 75, 174426 (2007). 

19 S. Mahmoudian and A. Langari, Phys. Rev. B. 77, 024420 (2008). 



20 H.Rezania, ALangari, P.Thalmeier, Phys. Rev. B 74, 094438 
(2008). 

21 F. Wegner, Ann. Phyzik. 3, 77 (1994). 

22 S. D. Glazek and K. G. Wilson,Phys. Rev. D 48 5863 (1993). 

23 S. D. Glazek and K. G. Wilson,Phys. Rev. D 49 5863 (1994). 



A. Mielke, Eur. Phys. J. B. 5, 605 (1998). 
Q. Gu, Phys. Rev. B 66, 052404 (2002). 



24 
25 

26 LA. Wu, M. S. Sarandy, and D. A. Lidar,Phys. Rev. Lett 93 
250404 (2004). 

27 G. Vidal, J. I. Latorre, E. Rico, and A. Kitaev, Phys. Rev. Lett. 90, 
227902 (2003). 

28 J. Vidal, G. Palacios, and R. Mosseri, Phys. Rev. A 69, 022107 
(2004). 

29 T. J. Osborne and M. A. Nielsen, Phys. Rev. A 66, 0321 10 (2002). 

30 I. Bose and E. Chattopadhyay, Phys. Rev. A 66, 062320 (2002). 

31 F. Verstraete, M. Popp, and J. I. Cirac, Phys. Rev. Lett. 92, 027901 
(2004). 

32 P. Zanardi and X. Wang, quant-ph/0201028 

33 Shi-Jian Gu, Shu-Sa Deng, You-Quan Li, and Hai-Qing Lin.Phys. 
Rev. Lett. 93,086402 (2004). 

34 Alberto Anfossi, Paolo Giorda, and Arianna Montorsi, Phys. Rev. 
B75, 165106(2007). 

35 T. R. de Oliveira, G. Rigolin, M. C. de Oliveira and E. Miranda, 



Phys. Rev. A. 77, 032325 (2008). 

36 H. Rezania, private communication. 

37 J. J. Mendoza- Arenas, R. Franco, and J. Silva- Valencia Phys. Rev. 
B 81, 035103 (2010); Phys. Rev. A 81, 062310 (2010). 

Appendix A: Effect of q n operators 



Appendix B: Spectrum of two and three dimensional model 

The excitation spectrum of the two dimensional model is 
given by the following equation where {k x , k y ) represent the 
momentum components, 



uj(k) = 1 + A( 2[cos(fc E ) + cos(k y )] + 8a ) — A ( — Arj z [cos(k x ) + cos(k y )] + cos(2/c K ) + 2[cos(fcj; + k y ) + cos(k x — k y )] 



cos(2fc, y ) + 2rj; + A - 



21rj 2 z + 23 



[cos(k x ) + cos(fcj,)] + 3[cos(k x + 2k y ) + cos(k x — 2k y ) + cos(2fc.£ + k y ) 



+ cos(2fc 2 , — k y )] + cos(3fc x ) + cos(3fc. y ) + (16a + 8rj z )[cos(k x + k y ) + cos(fc x — k y )] + 16rj z + 32a — 10^. 



(Bl) 
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TABLE IV. The q n effects on unperturbed bipartite eigenstates 
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The one magnon spectrum of the three dimensional model 
is given by the following expression: 

J 



w(fc) = l+A[2ai + 12a]-A 2 [37/ 2 2 -47/ 2 a 1 +2a3+a 2 ]+AV7 2 ^ ( B2 ) 

I 

which oi, a2 and 03 are Appendix C: Correlation functions 



ai = cos(fc x ) + cos(A; y ) + cos(fc 2 ), 
02 = cos(3fc x ) + cos(3fcj,) + cos(3fc 2 ), 
a 3 = cos(k x + k y ) + cos(fcj. — k y ) + cos(fc x + k z ), 
+ cos(k x — k z ) + cos(k y + k z ) + cos(k y — k z ). 

(B3) 



The correlation function of the nearest neighbor spins on 
the d-dimensional hypercubic KN lattice can be obtained via 
the generating functions defined in Eq.(l24l>. The correlations 
functions get a simple form for the isotropic case in the local 
interactions, a = (A = 1), which are presented in this ap- 
pendix. The correlation functions between the nearest neigh- 



12 



bor r spins are: 



H 'i+l 



+4(4d- l)r, x (r£ + nl + t?2) - I6(4d - 3)4], (CI) 



( T f T f+i) = - \vy X - g^%A 2 + ^[-8?/ y (^ + ri 2 z ) 
+4(4d - l)^(r,2 + v 2 y + V 2 Z ) - 16(4d - 3)^], (C2) 



A 
4' 



3A 2 



1G 



4?7 z x 2 - 
(n-1) 



^0* 



i 2 



i 2 ) 



* 2 te)(% 2 )], 



(C7) 



1 _ 3 _2 A 3 
(r?rf +1 ) = ——r) x \ - -VxVy* + ^[-^Vzivl + vl) _ _ 2 

A / 3a: A / / 

+4(4d - l)vz(vl + V 2 y + vl) - 16(4d - 3)^]. (C3) K^.iH -~n x —V V V Z 

3 

The analytic form for the local correlation functions be- A 2 



tween the impurity spin (s) and the spin on the lattice (r) are 



VM +%+ Vl) + 3x 4 r, x ( Vx 2 + Vy 2 + Vz 2 ) 



in the following form -16, 2 (^) _ 2x S x (v v 2 + _ 4*^], (C8) 



= "7 + l^v +? ?^ a2 ~ d^VxVvVz, ,vvs A ' 3a; 3 A 2 - - 



Id -2-3 ' l ' 4 

( T f s f> = -7 + t(^ 2 + ) A - dX VxVvVz, T 3 



^ 1 A Vlx < l lzJ" 'ix'iyizi \" 

+— [4~ 2 " 

= -7 + 7^* + ^ ~ ^VxVyVz- (C4) 



+-5-[4a^^(»S + ^ + »E) + + Vy + Vz) 



4 4 *' -16, 2 (^)-2A;(^ 2 +% 2 )-4xV], (C9) 

The correlation functions for the n-legged ladder are as the 
followings: 

1 x x \ A 3A 
{ T l,i T l,i+l) = -4% - —VyVz 

+ 7r[3^(^ 2 + ^ 2 + ?7 2 ) + 477.x 2 + ^ 2 + ^ 2 ) 

-2r^ 2 + ry 2 + 2r, 2 ) - l 6 (^a*(,, a ) (,£)], = - 

A 3 , , 

j 3 -f +Y^ x2 Vz(vl + V 2 y +vl) + 3x%(rix 2 + Vy +Vz 2 ) 

« i <i + i)= -4% - — -lto 8 ^) - 2^;(r,: 2 + r,' 2 ) - 4xM 3 ]. (CIO) 

—3 

+ ^[3%0? 2 + ^ + vl) + 4v y x 2 ^-^(Vx + Vy + Vz 2 ) 
~2r h ( v l + r,l + 2 V 2 ) - 16^—^-x 2 (r, y )(ri' 2 )}, (C6) 



